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SYNOPSIS 


ELECTRONIC STRUCTURE OF Cr-Al ALLOYS 


SURENDRA KUMAR VERMA 

Department of Physics, Institute of Basic Sciences, 
Bundelkhand University, Jhansi (India) 


The aim of this work is to study the electronic structure of disordered 
Cr-Al alloys. The electronic structure of this system cannot be studied by 
using standard method of band theory because disordered system lacks the 
translational invariance. During the past decades, Coherent Potential 
Approximation (CPA) theory has been very successful in understanding the 
electronic structure of disordered alloys. The CPA is a mean field 
approximation in which the disordered alloy is replaced by an effective 
medium, which is determined by a self-consistent condition. This theory 
when applied to muffin-tin potential of disordered alloy is known as 

R1-R5 

Korringa-Kohn-Rostoker Coherent Potential Approximation (KKR-CPA) 
Within the local density approximation of the density functional theory, it 
can be made fully charged self-consistent. Thus, the charged self-consistent 
KKR-CPA provides a first principles parameter-free theory of electronic 
structure calculation of disordered alloys. 






XI 


We have chosen Cr x Alj_ x alloys because it is important Body 
Centered Cubic (bcc) alloy system. The system forms solid solution for 

0.745 < X <\ in a bcc paramagnetic phase . To my knowledge, no such 
calculation has been reported for these systems till now, although some 

R7 

experimental results are available on Cr-Al alloys . We shall focus on the 
results of Density Of States (DOS), as this is the basic quantity needed to 
explain the various experimental results. 

We shall apply the charged self-consistent KKR-CPA method for 
calculation of DOS of Cr-Al system. Here, we shall outline some important 
equations involved in the method. Let us consider a binary alloy a x b y 

where X and Y are atomic concentrations. The basic quantity of interest is 
the ensemble-averaged Green’s function of the alloy. The Green’s function 

R1 

for a particular configuration of the alloy can be written as 

G(E,r,r')= (1) 

W L 

Here, r n and r m ' are within n ^ and m^ muffin-tin spheres and the wave 
functions Z”{£,rJ and are respectively the regular and irregular 

solutions of the Schrodinger equation for the single muffm-tin potential 
centered at r n , and the T nm are path operators. The ensemble-averages of 
equation (1) are related to the ensemble average of the path operators, which 
is determined by invoking the CPA. The CPA condition, that the average 
scattering from each site must be zero, can be expressed 

IfL = k-tf +!'■'? +(<? -<» L P) 
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Here t A(B) denotes on-the-energy-shell matrix element of the t-matrix of an 

isolated A(B) atom and C is used to label quantities for the CPA medium. 
The path operator matrix is given by 


h 0 , L = n/M ftfc'to- ^wrL^ 


( 3 ) 


Here, B R are the KKR structure constants and Q is the unit cell volume. For | 

full charged self-consistency, the charge densities in A and B cells are j 

needed. These can be obtained from the restricted site averages of the j 

Green’s function <G(E,r,r’)> A(B) as j 

. 

pA(B){r ) = -(l/^)£lm < G(E, r,r')> A(e} dE (4) j 

■ ■ ■ ■ ' ' ' . ■ . . . ■ . I 

where p A (r) and p B (r) denote the charge densities associated with A and B 

cells respectively. 

The component DOS for the A and B atoms are expressed as 

. : '' ' • ' ' X ' . ' : ■■ ' 

• ■ ■ ■ ' ■ ■ ; , ■ ■ ■ ■ : ■ ■ 

- . ' ' , ' : . ■ - ■. ' . : ‘.' 7 " ■- 7 ^ ■. ; ■■ ■ . : . : ■■ . . , ; ' ' 

I Pa(b)( e ) = -(0r)lm £< G(E,r,r') > A{B) dr (5) 

fj . ' 1 

The average density of states for the alloy is calculated as 
i p(E) = X.p A (E)+r-P,(E) (6) 

I 1 ' ■ . .. ■. ' ■■■■■ . ■■■ ' ■ ■ ■ ' ' ' ■ ■ ■■■ , 
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For the computational details of charged self-consistent KKR-CPA, 
we will start with atomic charge densities (p in A ,p in B ) for ,4 and B and will 
calculate the potentials {v‘\V‘") using the local density approximation. The 
^’(z)and t~\x) will then be computed for the input potentials ). We 

will use the tetrahedron method for the calculation of path operator 

After solving the KKR-CPA equation (2) iteratively, the new charge 
densities (pT>p7') will be obtained from equation (4). The new potentials 
(f“" , V,““ ) in the alloy will then be calculated. The whole calculation is 
repeated, until the difference of old and new charge densities is within a 
certain tolerance (integrated Root Mean Square (rms) difference between 
input and output charge densities (0.0001 electrons)). Fermi energy will be 
calculated by generalized Lloyd’s formula 118 , which does not give unphysical 
jumps in the integrated density of states. 

We shall calculate the DOS of Cr x Alj_ x for various concentrations 

(. X = 1 . 0 , 0 . 95 , 0.85 andQ.15) of Cr in the alloys. We shall compare our results 
of DOS with Cr-based rigid band model. We shall also study the DOS at 
Fermi energy (Ep) verses Cr concentration in alloys. We shall also compare 

R7 

our results of DOS with available experimental results . 

Since the KKR-CPA calculation takes a lot of computer time, 
therefore, we shall propose new scheme based on double peak semicircular 
model for the alloys. This proposed model would be modified form of earlier 

R9 

multi-peak model . In the earlier multi-peak model cotangent of phase 
shifts of either constituents was modeled with the help of step functions and 
it was assumed that their derivatives are zero at boundaries. But in the 







proposed scheme steps functions will not be used. We shall provide band- 
centers and bandwidths of constituents of alloys, and then we shall predict 
DOS of the alloy. The proposed model will also be the aim of the thesis. 
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INTRODUCTION: 

There are a rich variety of disordered alloys having diverse electrical, 
magnetic, transport, optical and super conducting properties in nature; these 
properties make them very useful for technological applications. For 
understanding these properties knowledge of their electronic structure is 
essential. In disordered alloys due to the lackness of translational invariance, 
which is a characteristic of ordered solids, Bloch’s theorem, which greatly 
simplifies the electronic structure calculation of ordered solids, is 
inapplicable to these systems. In other words, the standard band theory 
method developed for ordered solids cannot be applied to disordered alloys. 

The main subject of this thesis is the calculations of electronic 
structure of substitutional random binary alloys which are the simplest kind 
of disordered systems. In such alloys there exists a lattice structure but each 
lattice point can be occupied by either of the constituent atoms. Examples of 
such systems are a-brass Cu-Zn, Cu-Ge, Cu-Pd, Cu-Ni, Cu-Au, Ag-Pd, Nb- 
Mo, Cr-V and Cr-Mo etc. Such an alloy will be denoted by A X B Y , Where X 
and Y are concentrations of A and B types of atoms in the alloy, 
respectively. 

A great progress has been made in understanding the electronic 
structure of the disordered alloys during last three decades, by the 
application of Coherent Potential Approximation (CPA ) 1 ' 3 . The CPA is a 
mean field approximation, which for the sake of configurational averaging 
replaces a disordered alloy (A X B Y ) by an ordered solid of effective atoms. 
These effective atoms are determined by using the self-consistent condition 
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that if an A (B) atom is embedded in this effective medium, the average 
scattering with respect to the medium is zero i.e. 

x// A +y/ B M =0 (1.1) 

where 4™ is atomic scattering operator for an A or B atom embedded in the 

effective medium. A simpler approximation is average t-matrix 

approximation (ATA) 1 in which the t-matrix corresponding to the effective 
atom is 

t AU =x.t A + y.t B (1.2) 

where / A(B ) is atomic scattering operator for an isolated A(B) atom. Another 
simple approximation is virtual crystal approximation (VCA) 4 in which the 
potential corresponding to the effective atoms is assumed to be 

V VCA =x.V A +yV B (1.3) 


where Va(B) is the potential of A(B) atom. This approximation is good only 
when the difference between the two constituents potential is small. Note 
that all these approximations (CPA, ATA, VCA) are single-site 
approximations; i.e. they neglect correlated scattering from clusters of 
atoms. The CPA has been found to be the best single-site approximation for 
calculation of electronic structure of disordered alloys. 

A frequently quoted model in alloy theory is the rigid band model 5 . In 
this model the potentials of the constituent atoms are assumed to be the 
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same. However, Fermi energy (Ep) is adjusted to give the required number 
of valence electrons (N v ) per atom in the alloy as 


N v = \p{E)dE 


(1.4) 


where p(E ) is the density of states (DOS) per atom of the host system. The 
rigid band model is a crude approximation and does not explain many 
experimental results. Because of its simplicity, it is frequently used as a first 
step to get an idea of the electronic structure of the alloy. 

There is very large number of interacting electrons (« 10 23 ) interacting 
with nuclear potentials in a solid. Therefore, the calculation of energy levels 
of such a system is a many-body problem, which cannot be solved exactly. 
However, this problem can be solved to a great degree of accuracy by using 
the one electron approximation. In this approximation, the many-body 
problem is reduced to a problem of one electron moving in some effective 
potential which is determined self-consistently. Muffin-tin approximation is 
further used to simplify the calculation of this potential. In this 
approximation the potential is assumed to be spherically symmetric within a 
sphere and Wigner-Seitz cell as 


V{r). 


V v) (r) 


for r < r 


yl 

y A(B) 


for r > r 


(1.5) 


The V 1 is a constant interstitial potential as 6 


v 1 

V A(B) 


3 }v{r)r 2 dr Y (r v 


r 3 ) 

' m / 


( 1 . 6 ) 
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where V(r) is a spherically symmetric potential, which, it is hoped, will be 
slowly varying in the region between the muffin-tin sphere and the Wigner- 
Seitz sphere of radius r m . For a pure A(B) solid. Vj (B) defines muffin-tin 

zero, i.e. this constant is subtracted from the potential (1.5) making 
interstitial potential zero. For the alloy, this is calculated as 

Kuor = x-vl + y-v» ( 1 . 7 ) 

Note that muffin-tin potentials are non-overlapping. The muffin-tin 
approximation is a reasonably good approximation for the potential as can 
be seen by referring to the book by Moruzzi et al 7 . Using this approximation, 
they have calculated various properties such as cohesive energy, bulk 
modulus, density of states and Fermi energy. They got very good agreements 
with experimental results and were able to explain general trends. Recent 
experience with alloys show that the muffin-tin approximation is also a good 
approximation for metallic alloys 1 ' 3, 8 ' 14 . The muffin-tin potentials for the 
constituent atoms are constructed by using the local density approximation 
(LDA) of the density functional theory 15 . The potential can be written as 3 


Vm (?) = V3B. (f) + V N W (?) + F* S) (?) (1.8) 

where coulomb, nuclear and exchange-correlation contributions are 

P ) 


\df 


r - r 


( 1 . 9 ) 



( 1 . 10 ) 
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( 1 . 11 ) 


V A(H) (r) = ^ x< \-P- A{R) ^\ 

XC Sp A(B) {r ) 

where p A(H) is the electron charge density in A(B) cell, Z is the atomic 

number and Q denotes the integral over the Wigner-Seitz unit cell. Here, 
Exc is the exchange-correlation energy functional, which in the local density 
approximation is given by 15 

E X c {p} * \p(r)e xc (i p(f))dr (1-12) 


where e xc {p) is the contribution of exchange and correlation to the total 
energy (per electron) in a homogeneous but interacting electron gas of 
density p . For paramagnetic case, von Barth-Hedin form of s xc is 16 


s B xc (r s ) = - 0,91633 -0.045f( 


21 


(1.13) 


where r s is given by 


4 ^, ,3 


W = 


f-1 

Ip, 


and function 


F{x) = (l + x 3 )/«(l + l/x)-x 2 +x/2-l/3 


(1.14) 


(1.15) 


Gunnarsson and Lundqvist form of s xc for paramagnetic case is 
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£< xc( r s ) : 


0.458 


r v \ 


-0.0333F 


V 11.4y 


(1.16) 


Vosko etal. form of s xc for paramagnetic case is 
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£ xc ( r s > : 


0.91633 
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' x 2 ' 

2b , 

_l fan * 

r q i 


_ X(x)_ 

i tall 

Q 

_2 x + b _ 


bx n 


X(x 0 ) 


In\ 


( x ~x of 

X(x) 


+ ^ + 2x o) tan -i 
Q 


Q 


2 x + b 


where 


X(x) = x 2 +bx + c , Q = (4c-b 2 ) 1/2 and x=(r s ) 1/2 (1.17) 

Here, A = 0.0621814, x 0 = -0.10498, b = 3.72744 and c= 12.9352 

It has been seen that use of different exchange-correlation does not 
give large differences in the calculation of various properties 15 . 
For example differences in cohesive energies in Li, Na, K and Rb were 
found to be less than 8 mRy 18 . Because Moruzzi et al 7 , have found good 
agreement between theory and experiment for several metals including Cr 
and Al, using von Barth-Hedin form of the exchange-correlation potential. 
Therefore, We have used this form in our calculations for above reason as 
will be discussed in detail in Chapter 3. 

The CPA theory has been very successful as a single-site 
approximation to calculate the electronic structure of random substitutional 
disordered alloys during last thirty years. This theory when applied to 
muffin-tin model of disordered alloys, is known as Korringa-Kohn- 
Rostoker-CPA (KKR-CPA) theoiy 1 . In the perfect crystal limit, the 
KKR-CPA theoiy reduced to the standard KKR band theory 1 . In earlier 
stage of application, potentials V A and V B were not determined self- 
consistently. In the most sophisticated application of the KKR-CPA, within 
the framework of the local density approximation V A and V B are determined 
self-consistently described in Chapter 2. This fully charge-self-consistent 
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KKR-CPA theory is a first principle parameter-free theory of the electronic 
structure of random alloys. 

The outline of the thesis is as follows. In Chapter 2, we will present 
the brief formulation of charge self-consistent KKR-CPA. Section 2.1 
presents the introduction of the chapter. In section 2.2, we derive the KKR- 
CPA equation in multiple scattering formalism. In section 2.3, we present 
the charge self-consistent KKR-CPA method. 

In Chapter 3, we will apply the charge self-consistent KKR-CPA 
method in calculation of electronic structure of Cr-Al alloy. Section 3.1 
present the introduction of the chapter, we will give the computational 
details of our calculation in section 3.2. The Green’s function method is used 
to calculate charge densities of constituent atoms in the alloy. The flow chart 
of charge self-consistent KKR-CPA is also given in this section. The mixing 
scheme for CPA and charge self-consistency loop are discussed. The study 
of charge self-consistent potentials, charge densities and the density of states 
as a function of Cr concentration in the Cr x Al Y for X= 0.95,0.85 and 0.75 
are given in section 3.3. In this section we also compare the results of charge 
self-consistent KKR-CPA density of states with the density of states using 
Cr-based rigid band model. We have chosen von Barth-Hedin exchange 
correlation potential to calculate the new crystal potentials of constituent 
atoms in the alloy. For Cr-Al alloy, we compare the results of density of 
states calculated by using the charge self-consistent KKR-CPA method with 
the X-ray photo spectra experimental results. We find that there is a very 
good agreement between our results and experimental results for different 
concentration of alloy. It is also found that Cr-based rigid band model fails 
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in the calculation of electronic structure of Cr-Al alloys. In the final section 
3.4 we summarize our main conclusions. 

In Chapter 4, we apply the KKR-CPA to the s-phase shift double peak 
semicircular model. Section 4.1 provides an introduction of this chapter. In 
section 4.2 we discuss the s-phase shift semicircular model in the KKR 
framework. In section 4.3, the path operators and t-matrices for the s-phase 
shift double peak semicircular model is derived. In section 4.4, formula for 
the KKR-CPA density of states to the s-phase shift double peak semicircular 
model is derived. In section 4.5, we present the CPA equation and 
computational details for this model with programme. In section 4.6, we 
discuss our results for s-phase shift double peak semicircular model. In 
section 4.7, we summarize our main conclusions. 

Finally, in Chapter 5, we give the summary of the achievements of the 
present work. Some suggestions for further improvements are also given and 
some possible remaining problems are mentioned. 


CHAPTER - 2 


CHARGE SELF-CONSISTENT KKR-CPA 


2.1 INTRODUCTION 


To understanding the electronic structure of disordered alloys, the 
Korringa-Kohn-Rostoker Coherent Potential Approximation (KKR-CPA) 
theory has been very successful in last three decades as was mentioned in 
chapter- 1. It is basically an application of the coherent potential 
approximation (CPA) 1 ' 3 to muffin-tin potential of the disordered alloys. In 
the CPA, one replaces a disordered alloy by an ordered array of effective 
atoms which are determined by a self-consistent condition. The muffin-tin 
approximation of the potential plays an important role in simplifying the 
theory. Because the constituent muffin-tin potentials are spherically 
symmetric within a certain radius r m , various quantities such as t-matrices 
can be expressed in simple forms in the angular momentum space. The non- 
overlapping nature of the constituent potentials further simplifies the 
multiple scattering equations to tractable forms. 

The KKR-CPA reduced to standard KKR band theory in the perfect 
crystal limit. Within the local density approximation of the density- 
functional theory (DFT) 15 " 18 , it can be made fully charge self-consistent. 
Thus charge self-consistent KKR-CPA provides a first principles parameter- 
free theory of electronic structure of disordered alloys. The theory rests on 
the same theoretical footing as the band theory of pure metals. 

The outline of this chapter is as follows. In section 2.2, the 
formulation of the KKR-CPA will be briefly discussed. In section 2.3, we 
present a brief formulation of the charge self-consistent KKR-CPA. The 
final section 2.4 gives conclusions. 
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2.2 KORRINGA-KOHN-ROSTOKER COHERENT-POTENTIAL- 
APPROXIMATION (KKR-CPA) 

Consider a disordered substitutional binary alloy A X B Y of A and B 
atoms with concentrations X and Y respectively. The one electron potential 
in the alloy can be written as 1 " 3 


m = ZK(n) ( 2 . 1 ) 

/=l 

where the vectors r i are defined by r t =r- R t ,R, are the location of the lattice 
points and N is the total number of sites in the solid. V a (F ) (a=A,B) is the 

potential due to A or B atom and is assumed to be the muffin-tin form and 
does not overlap with each other; 

K(n) =v a(ri) , f < r m 

V a ( 7 , ) = Constant, rj > r m (2.2) 

where r m is the muffin-tin radius. 

The Hamiltonian for a certain configuration of the alloy in atomic mass unit 
is 

H = -V 2 +V(r) (2.3) 

The Green’s function in operator notation for a system of scatterers is 
defined as 

G = (EI- H)~ l (2.4) 
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where E and I are energy and identity matrix. The equation (2.4) can be 
expanded using the Dyson’s equation 1 as 

G = G 0 +G 0 TG 0 (2.5) 

where G 0 is the free electron Green’s function and 

T = F + FG 0 T (2.6) 

is the total scattering operator. T may be written as 

T = IX (2.7) 

iJ 

where the operators T t J are called path operators and satisfy the following 
multiple scattering equation 


Tij - ij + t,G 0 ^ T kj 

k*i 

Here, 

/,=F,(H-G 0 /,) 


( 2 . 8 ) 


(2.9) 


is the t-matrix that describes the scattering from an isolated potential on the 
i-th site. Equation (2.9) in angular momentum space has on-the-energy-shell 
matrix elements as 

= ~W l expO'^)sin(^) (2.10) 

where / is the angular momentum index, z = ( E ) U2 an d 5, are phase shifts. 
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The free electron Green’s function Go in real space can be expressed as 2 

L 

LL ' 

( 2 . 11 ) 

where r=r n +R n , r'= f m +R m and r < (r > ) is the smaller (greater) of the 
variable r n and r’ m . Also L= (l,m) is a composite index, Y L (x) is the real 
spherical harmonic and J L and h L are the spherical Bessel and Hankel 
functions respectively. The matrix B is defined as 


BZ’iz) = 


4«. zE y,, A-*.)* K (z 




- l\ 

; 


o-o 


where 


Cu. t| = 


( 2 . 12 ) 

(2.13) 


and h] denotes the outgoing Hankel function. 

Starting from equation (2.8) and after some algebraic manipulations, on the 
energy-shell matrix elements of the path operator can finally be expressed as 

T“'(Z)=[k'(T)-B(Z)F 1 |l (2.14) 

The Fourier transform of equation (2.12) 
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(2.15) 


k«L = 0/ aos = xp [- » A - a.)ki'(z) 

nm 

is related to the well-known KKR structure functions 19 . Equations (2.10) and 
(2.15) are very important in the KKR theory; t(%) depends only on the 
potential while B^(%) depends only on the lattice structure. These two 

matrices play a central role in the KKR-CPA theory. 

Now the Green’s function G can be written as 2 



(2.16) 

where r and r ' are within n th and m th muffm-tin spheres and the wave 
functions Z" im) (E,r n{m) ) and Jf m \E,r n{m) ) are respectively the regular and 

irregular solutions of the differential equation. 


[ -^ 2 +K,,,(r nW )-E ]z; -0 


(2.17a) 


For spherically symmetric potentials one can write 


Zf w (E,r) = Y L (r)Z?*\E,r ) 


(2.17b) 


Jf} B) {E,r) = Y L (r)af B) (E,r) J 

where Zf {B) and af w are radial wave functions and are normalized such that 
for r > r m they join smoothly to 
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Z L*\z>r) ~ J iXx,r)[t'; w ]' -i.%.h L (x,r) 

oc MB Xx,r) = J L { X ,r) 


(2.17c) 


The regular wave function Z AW can be written as 20 

Zf w =^ B XE)ix A L iB) (r,E) (2.17d) 

where <fi A{H) (E ) is an energy dependent renormalization factor independent of 
r such that y/ A{li) -» r 7 ' for r -» O 20 . 

Equations (2.16) and (2.14) are exact within the muffin-tin 
approximation and can be used to obtain the Green’s function for a cluster of 
atoms 21 . However, here we are interested in an infinite system and would 
like to obtain ensemble average of G(E,r,F ) . The ensemble averages of 
equation (2.16) are related to the ensemble averages of the path operators, 
which will be determined by invoking the CPA. The CPA condition that the 
average scattering from each site must be zero, can be expressed in terms of 
Green’s function operator as 

x. < G > A +y. < G > B = G (2. 1 8) 

where < G > A{H) denotes the restricted site average of G when zeroth site is 
occupied by an A (B) atom. Equivalently this condition can be written in 
terms of path operators as 

*• < T nn >n*A +y- < T m >«=B= T m ( 2 - 19) 
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where < T nn > n=A and < T m > n=R are the restricted site averages of T nn where 
n th site is occupied by an A (B) atom and T n c n is the nn path operator for the 

CPA medium. After solving equation (2.19), we get the KKR-CPA 
condition as 

4A.L + )'[A,L = 'u- (2.20) 

where 

K«L=[ ;+r » ««>-'?)£ (2.2i) 

and C is used to label quantities for the CPA medium. The path operator 
matrix is giver by 

fcl, = (2.22) 

where Q is the unit cell volume. We note that involves a complicated 
Brillouin zone integration. Though B R (%) is an off-diagonal matrix, for cubic 

symmetry the integral reduces to a diagonal matrix 3 for 1 < 2 . By simple 
manipulations, one can reduce the KKR-CPA equation (2.20) with the help 
of equation (2.21) to a computationally simpler form as 

Iff; 1,/; = \. X ^A +y*B + 0c ~*A )T 00 Oc ~*B )]u' (2.23) 

This equation can be solved only by a numerical iterative method. The input 
for iterating this equation is r A \, t~ B ] via equation (2.10) and B R (%) via 

equation (2.15). 
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2.3 CHARGE SELF-CONSISTENT KKR-CPA 

For full charge self-consistency, the charge densities in A and B cells 
are needed. These can be obtained 2 , once the restricted site averages of the 
Green’s function < G(E,r,F ) > A(B) are known as 

Ef 

Phb)(?) = ~(/y Jim <G(E,r,F)> A{B) dE (2.24) 


where p AW denotes the charge density associated with an A (B) cell. 
Equation (2.24) can be further rewritten as 

- -( :%) flm tr[z m (r)Z m (r')D m) 'E\ LL .dE (2.25) 

-00 

For the charge self-consistent KKR-CPA, the KKR-CPA equation (2.23) is 
first solved for given potentials of A and B atoms. The new charge densities 
p A and p tt are calculated from equation (2.25). New potentials are then 
calculated by using the local density approximation of the density-functional 
theory as has been described in Chapter 1 . The new potentials are used to 
solve the KKR-CPA equation and again a new set of potentials is calculated. 
We iterate this process until the potentials and charge densities get 
converged. 

To calculate the density of state and component density of states, we 
define the matrices F A and F B as 2 

FfP = \zr (r)Zt- B) (rW (2.26) 

ci 
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where Cl denotes the integral over the unit cell. Then the component density 
of states for the A and B atoms can be expressed as 


Pa W W = ~(%) Im J< G(E,r,F) > A(B) dr 

Cl 

or 

= -(j/)Im tr[F Am D im T^\ L . (2.27) 

where tr denotes the trace in L-space. Since equation (2.27) is diagonal in 
L-space, it allows L decomposed symmetry components (s,p,t 2g and e g ) of 
Pa(B)(E) i n the alloy. The cubic symmetry reduces the number of distinct 

elements to four, i.e., s,p,t 2g and e g . The average density of states for the 
alloy is calculated as 

P(E) = +y.F'B,teL, (2.28) 


2.4 CONCLUSIONS: 

In this chapter basically we have shown main equations of charge 
self-consistent KKR-CPA method. The detailed computational procedure 
and the flow chart will be given in Chapter 3. 
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1 CHAPTER- 3 

:: :: 'y£x: ; - v - ’* , : ; V;^ . 

Is a APPLICATION OF CHARGE SELF- CONSISTANT 


KKR-CPA TO Cr-AI ALLOYS 


3.1 INTRODUCTION 

In this chapter we will calculate the electronic structure of disordered 

paramagnetic Cr-Al alloys by the application of the charge self-consistent 
'22 

KKR-CPA . We have chosen this system because it is very important bcc 
alloy system. Cr x Al Y alloy system form bcc solid solution for 1 > A > 74.5 23 . 
To our knowledge, no such calculations have been reported for this system 
till now, although some important experimental results of X-Ray photo 
spectra (XPS) and X-Ray emission spectra (XES) of the valence band on 
concentration X=1.0, 0.95, 0.85 and 0.75 are available on Cr-Al alloys 24 . In 
experimental results, the shape and half width of XPS of valence band 
changes notably at X=0.85 and reported largest at this concentration. In this 
chapter we shall focus on the results for the density of states, as this is the 
basic quantity needed to explain various experimental results. We shall 
compare our calculated results of DOS with XPS of valence band of Cr x Al Y 
at X=0.95, 0.85 and 0.75. We shall calculate the half width of DOS band and 
show its theoretical behavior. 

The outline of this chapter is as follows. In section 3.2 we give 
computational details involved in the method with flow chart. In section 3.3, 
we present the results of electronic structure of Cr-Al alloys by using the 
charge self-consistent KKR-CPA method. We present our results of density 
of states as a function of Cr concentration in Cr-Al alloy. We have compared 
charge self-consistent KKR-CPA density of states (DOS) results with the 
rigid band model density of states results. We have compared the density of 
states with the available experimental data for Cr-Al alloy and have found 


18 


very good agreement between the two. The final section 3.4 gives 
conclusions. 

3.2 COMPUTATIONAL DETAILS 

We have applied the charge self-consistent KKR-CPA method to the 
Cr x Al Y alloy for different concentrations of Cr (x = 1 .0,0.95, 0.85<m/0.75) . This 
alloy has bcc structure and their lattice constants are presented 23 in table 3.1. 


TABLE 3.1 


Lattice constants of Cr x Al Y alloys 


S.No. 

Cr concentration (X) 

■ 

Lattice constants (a.u.) 

1 

1.00 

5.453 

2 

0.95 

5.484 

3 

0.85 

5.508 

4 

0.75 

5.553 


For the charge self-consistent KKR-CPA, we start with atomic charge 
densities (p",p%) for Cr and Al, and calculate the potentials (Fj,F;")from 
equation (1.8) as explained in Chapter 1. The t A (x) and t~ B \x) are then 
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computed for the input potentials (V‘ n ,V‘ n ). The KKR-CPA equation is 
solved iteratively within a tolerance of 0.0001, to calculate the . Note that 
in each KKR-CPA iteration T ^ has to be calculated using equation (2.22), 
which involves a complicated Brillouin Zone integration. In the initial stages 
of implementation of the KKR-CPA theory, this was a major problem, 
which was eventually solved by using special direction tachnique 25 . 
Although this technique is quite simple, it does not give as accurate results 
as tetrahedron technique which has been developed by Kaprzyk et al 26 , for 
disordered alloys. Since, higher accuracy is essential for full charge self- 
consistency, we have used the tetrahedron method using 858 K points in the 
irreducible part of the Brillouin Zone. We have observed that it is difficult to 
achieve convergence if we use less number of K points in the tetrahedron. 
The new charge densities ( p , p ° B ut ) are obtained by equation (2.25). 
The new potentials (V™‘ , V ° ut ) in the alloy are then calculated using equation 
(1.8) of Chapter 1. Again, these new potentials are taken as input to repeat 
the whole process, until the difference of old and new charge densities is 
within a certain tolerance (e CHARGE ) . This tolerance has been defined as 

r m 

P ) __ f ( ry out — q * n ^ y ^ r!y f 3 , Ilf 

S CHARGE ~ J \Pa(B) r A (B) ) r ur ' 1 ) 

0 

This tolerance is taken to be 0.0001. The flow chart of the charge self- 
consistent KKR-CPA has been given below. 


20 


FLOW CHART OF CHARGE SELF-CONSISTENT 

Input: atomic numbers (Z A ,Z B ), lattice constant a, lattice structure and 
concentration (X) of A in A X B Y 


Guess charge densities (p A ,p%) 


Calculate input potentials (V A ,V‘ n ) 




Oc l ) ,n = x.r ; 1 + y.t-y 

^ / 


Obtain by KKR-CPA equation 


Is -(/- 1 )"’ < o.oooi 



Calculate E F , p™[ H) and V A 


Calculate p (E),p A (E) and p B (E) 



Chart-3.1 

2 


Qc) =(tc)° 




















For the new charge densities, we have to calculate Fermi energy, as 
can be seen from equation (2.25). The Fermi energy is calculated using 
Lloyd formula for integrated density of states. 

The Lloyd formula 1 is written as 


N , (E) = NUE) + -lm\ 

K 


x.lndet 



+ y. IndeU" 1 



1 _ 

izN 


Im^lndetj^' -B R 

K 


(3.2) 


where N'(E) is the integrated density of states and N‘ 0 (E) is the free electron 
integrated density of states. Though this formula is easy to handle, it 
sometimes gives unphysical jumps, which can occur in In (logarithm) of a 
complex quantity, used in the formula. Kaprzyk and Bansil 20 have derived a 
generalized Lloyd formula, which does not give such unphysical jump and is 
given below: 


N‘ {E) = (lml7t)tr 


! 7 £ln[G„-'(£,£) + J 0' -pc]' -ln(G c ) 


where 


+ x.ln[(^ (r m )) _1 ]+ y.ln[^ B (r m ))'‘ G g ] 


(3.3) 




(3.4) 


Here, the KKR-CPA Green’s function G c and the associated log-derivative 
J3 C are the self-consistent solutions of the CPA equation (2.18) 

x.G A +y.G B = G c (3.5) 
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where 


G < =^Z[G.-'(£^) + /8 y -yS‘ : ]■' (3.6) 

P? w {E) = r 2 ~\nZ A L (B \E,r) (3.7) 

and p[ is the log-derivative for the spherical Bessel function. Note that 
^ w (r)and Zi {H) (r) have been already defined in equation (2.17). 

Although equation (3.3) also contains the logarithm functions, but 
unphysical jumps of these functions cancel with each other. The integrated 
density of state at the Fermi energy gives the total number of valence 
electron in the alloy. Therefore, the Lloyd formula permits an evaluation of 
the Fermi energy without requiring an explicit computation of the density of 
states. We have used formula of Kaprzyk and Bansil, equation (3.3), to 
calculate the Fermi energy. 

The calculation of charge densities from equation (2.25) using 
integration along real energy is very slow and time consuming. To speed up 
the calculations, the complex energy method 27,28 has been used to calculate 
the charge densities and potentials. The K -space integration is much faster 
in the complex plane because the integrand becomes smoother when the 
energy is complex. 

To obtain full charge self-consistency it takes a large number of 
iterations to achieve convergence. Therefore, to facilitate both the CPA and 
charge self-consistency loops, we have used mixing schemes 7,29 . In the CPA 
loop new is taken as 

[(h' 1 ) W j/7?7M+i = a CFA l(L ) \rrcPA + (1 ~ a c:pA )l(fc ) hrcPA (^.8) 
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where a CPA is the mixing parameter for the CPA and ITCPA is the number 
of iterations. In the charge self-consistency loop, the new potential is taken 
as 

Vm)\,TCHAWH + \ = a MB)\FA%)\ rrCHAR GI. + 0 " a A(B))\F i,B) J ITCHARGE ^ 

where a v A(B) is the mixing parameter for A(B) potential and ITCHARGE is 
the number of iteration in charge density loop. We have taken a A(B) between 
0.75to 1 .00 for all concentrations of Cr in the alloy. The CPA mixing a A[B) is 
taken as 0.05. 

The CPU time on Silicon-Graphics Indy-R4400 system for running 
single charge self-consistency loop is about 6 minutes. The required number 
of iterations for the convergence of the charge densities in the charge self- 
consistency loop is about 400 for all concentrations of Cr in the Cr-Al alloy. 
The final calculation of density of states and component density of states is 
done along the real energy axis and takes about another 15 minutes of CPU 
time on this machine. 

Table 3.2 gives the Fermi energies (E F ) for various Cr-Al alloy 
composition and alloy density of states at the Fermi energy calculated within 
the charge self-consistent KKR-CPA. 
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Table- 3.2 


Density of states (DOS) at Fermi energy (E F ) for Cr x Al Y alloys 


DOS at E F (St/Ry.atom) 


0.75992 


0.74931 


0.74194 







3.3 RESULTS AND DISCUSSION 


We have applied the charge self-consistent KKR-CPA method to 
calculate the electronic structure of paramagnetic bcc Cr x Al Y alloys. We 
have calculated the charge self-consistent potential, charge densities, 
component density of states and density of states of Cr x Al Y alloys for 
concentrations X=0.95, 0.85 and 0.75. We have presented component 
density of states, density of states and density of states at Fermi energy for 
above concentrations, in figure 3.1 to 3.10. These results can be used 
directly for calculating electronic structure (Fermi surface, band structure, 
conductivity, optical properties, transport properties etc.) of the alloys. We 
have also compared our calculated density of states with experimental 
measured X-Ray photo spectra (XPS) of Cr x Al Y alloys 24 at concentration 
X=0.95, 0.85 and 0.75. 

Figure 3.1 shows the comparison of the average density of states 
calculated using the charge self-consistent KKR-CPA, with those of Cr- 
based rigid band model for X=0.95. We note that, although it is Cr-richer 
alloy, the Cr-based rigid band model does not work very well in this case. 
We notice that, there is no qualitative and quantitative difference between 
the two results. However there is a shift in the charge self-consistent KKR- 
CPA density of states towards the lower energy region with respect to Cr- 
based rigid band model density of states. This shifts of density of states is 
more at higher energies. 
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ENERGY ( RYDBERG) 


Figure. 3.1. The average density of states by using 
charge self-consistent KKR-CPA method (solid curve) 
and density of states obtained by Cr-based rigid band 

model (broken curve) for Cr 095 A / 0 05 . 


I 




Figure 3.2 shows the comparison of the average density of states 
calculated using the charge self-consistent KKR-CPA, with those of Cr- 
based rigid band model for X=0.85. We observed that, although it is Cr-rich 
alloy, the Cr-based rigid band model does not work in this case. We notice 
that, there is a marked qualitative and quantitative difference between the 
two results. 


Figure 3.3 shows the comparison of the average density of states 
calculated using the charge self-consistent KKR-CPA, with those of 
Cr-based rigid band model for concentration X=0.75. There is big qualitative 
and quantitative difference between the two results. We find that it is also 
Cr-rich alloy; but the Cr-based rigid band model totally fails in this case. On 
the basis of above results it is found that Cr-based rigid band model fails in 
the calculations of electronic structure of Cr-Al alloys. 


In figures 3.4, 3.5 and 3.6, we present the component density of states 
of Cr and A1 in Cr x Al Y alloys for concentration X=0.95, 0.85 and 0.75 
respectively. One can use these results in the analysis and calculations of 
electronic structure of these alloys. The component density of states can be 
compared with experimental results of soft X-Ray (emission and absorption) 
spectra. Therefore, these results will be very important for experimentalist. 
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iti- * 1 


DENSITY OF STATES ( STATES/RYDBERG-ATOM) 



Figure 3.2, The average density of states by using 
charge self-consistent KKR-CPA method (solid curve) 
and density of states obtained by Cr-based rigid band 

model (broken curve) for Cr U5 Ai 0A5 . 
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Figure 3.3. The average density of states by using 
charge seif-consistent KKR-CPA method (solid curve) 
and density of states obtained by Cr-based rigid band 

model (broken curve) for Cr 075 ^/ 025 . 
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Figure 3-4- The component density of states of Cr (solid 

curve) and Al (broken curve) in Cr aK Ai U5 . 
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Figure 3.5. The component density of states of Cr (solid 
curve) and Al (broken curve) in Cr 0 M Ai U5 . 




m 


3501 


CrinCr 75 Al 25 
Al in Cr.-75Al.25 




. ] 

iBil 

Illi 

H 

I i ‘iM 

1 

: 1 


■I 


: 




o 0 
Q t— 

> < 


Q § 

H- IT 
Z ^ 

lu m 

=7 UJ 

O t 

OL < 


300- 


250- 


200 - 


150 


100- 




0-8 -0-6 -0-4 -0-2 00 0*2 

ENERGY (RYDBERG) 


1 


■I 


H 


Figure 3.6. The component density of states of Cr (solid 
curve) and Al (broken curve) in Cr 01 ^Ai 0 25 . 
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In figure 3.7, we show the comparison of the charge self-consistent 
KKR-CPA density of states (our theoretical results) with available results of 
X-Ray photo spectra (XPS) 24 for Cr 095 Al 005 alloy. The structures (the position 
of peaks) in our theoretical results are well reflected in X-Ray photo spectra 
(XPS) results. The theoretical half-width of density of states (DOS) curve 
below Fermi energy is 2.20 Rydberg approximately. We find that our results 
are in good agreement with experimental results. 


In figure 3.8, we show the comparison of the charge self-consistent 
KKR-CPA density of states (our theoretical results) with available results of 
X-Ray photo spectra (XPS) 24 for Cr OS5 Al OA5 alloy. The structures in our 

theoretical results are very well reflected in X-Ray photo spectra (XPS) 
results. The theoretical half-width of density of states (DOS) curve below 
Fermi energy is 2.3 1 Rydberg approximately. It is noticed that the half-width 
of the band increases about 5% of the O 095 T/ 005 alloy, which is quite 
satisfactory with XPS experimental results. Thus we find that our theoretical 
results are in very good agreement with experimental results. 
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alloy. Thus we found that half width of the DOS band changes 
approximately 5% at 85 atomic percentage of Cr. It is less for 95 and 75 
atomic percentage of Cr, as it was found in XPS results 24 . Therefore, our 
results of calculated KKR-CPA density of states (DOS) are in very good 
agreement with XPS experimental results. 
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Figure 3.7. The average density or surw uy u=n. a 
charge self-consistent KKR-CPA. method (solid curve) 
and experimental X-ray photo spectra results (broken 

curve) for Cr m Ai m . 
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Figure 3.8. The average density of states by using 
charge self-consistent KKR-CPA method (solid curve) 
and experimental X-ray photo spectra results (broken 

curve) for Cr OS5 A/ OA5 . 
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Figure 3.9. The average density of states by using 
charge self-consistent KKR-CPA method (solid curve) 
and experimental X-ray photo spectra results (broken 
curve) for Cr 015 Ai 025 . 
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3.4 CONCLUSSIONS: 


We have applied the charge self-consistent KKR-CPA method to 
CrxAly disordered alloy system. We have calculated the self-consistent 
charge densities, self-consistent potentials, component density of states and 
density of states for concentrations X=0.95, 0.85 and 0.75 by using the 
charge self-consistent KKR-CPA method. It is found that charge-self 
consistent KKR-CPA density of states for different concentrations of alloy is 
in very good agreement with reflected X-ray photo spectra (XPS) results. It 
is found that half width of DOS band of these alloys increases approximately 
5% at 85 atomic percentages of Cr with respect to 95 and 75 percent of Cr, 
which explains the change in shape and half width of experimental XPS 
results. It is also found that Cr-based rigid band model fails in the 
calculation of electronic structure of Cr-Al alloys. 


ft, : 










CHAPTER- 4 

APPLICATION OF KORRINGA-KOHN-ROSTOKER 
COHERENT-POTANTIAL- APPROXIMATION TO 
S-PHASE SHIFT DOUBLE- PEAK SEMICIRCULAR 

MODEL 




4.1 INTRODUCTION 




Because of the lengthy computation and complexity involve in the 
KKR-CPA method, it is difficult and time consuming to calculate the 
electronic structure of realistic systems. The density of states (DOS) is the 
basic quantity to understand the electronic structure of solids. The DOS of 
realistic system have many peaks (structures) when we plot it as a function 
of energy. We have proposed the double-peak semicircular model by which 
one can estimate about the realistic density of states without expending more 
time and computational efforts. The S-phase shift semicircular model has 
several advantages, (i) Because of semicircular modeling, the involved k- 
space integration, required to obtain the site-diagonal path operator is 
bypassed, thus reduces the computational efforts, (ii) The KKR-CPA 
calculation for single peak model has been tried out before Soven et al. , 
Rajput et al . 31 and Yadav 32 . In the S-phase shift model we shall assume that 
the S-phase shift dominates and B, t and T all the matrices involve in the 
KKR-CPA method reduces in the scalar form in the angular momentum 
space. The CPA equation for this model can be solved transparently by 
reducing it in polynomial form and using the computational method. 
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4.2 THE S-PHASE SHIFT SEMICIRCULAR MODEL 

For this model we have assumed that the phase shift has resonance 
behaviour and the cotangent of phase shift C A (b) is as 


where Ea(B) is the resonant energy (band center) and T a(B) is half- width of 
resonance. Another approximation is that k-space structure function is 
assumed as 


has no explicit dependence of energy (i.e. independent of energy) where 
% = E u 2 and B- q is the Fourier transformation of real space structure function 

B nm . The B- q { X ) is as 


In equation (4.3) by translational property one index m is fixed and we sum 
only on index n, then on the shell energy t-matrix element 
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<X) = -(%)-' e‘\ Sin(S L ) 


(4.4) 




where S L is the phase shift. 

From equation (4.4) we can write 


Z-CotS + i.% 


If we write Cot{5 ) = C , then 


Now the explicit scalar form of path operator for the model has been 
derived as follows 


with the help of equation (4.2) 


We can write 
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If we take the distribution of B- as 


From above equation we can see 


We will choose the distribution function f°(b) which must satisfy the above 
condition thus it can be written as 


When \b\ < 1 


f\b)db 


When \b\ > 1 


f\b)db 


X.\ *(C + ») 

By substituting b = Cos 6 it will reduces as 



n.% 0 J (C + Cos 0) 


k.x J 2Z\Z 2 +2CZ + 1) 


where /(Z) 


2Z 2 (Z 2 +2CZ + 1) 


The poles of integrands are Z =0 of order 2, other two poles Z=Zi and Z=Z 2 
are the roots of equation Z 2 + 2CZ +1=0, such that 


We can find Z x and Z 2 by solving equation (4.14) and (4.15) as 



(a) The Residue at Z=0 of /(Z) 


Limit 


(b) The Residue at Z=Z\ and Z=Z 2) if we take + is imaginary part in the 
energy then 


E = E + is 


which will be as 


for 8 -* 0 


Now from equation (4.14) if |Z,| > 1 then |Z 2 | < 1 and hence Z, will lie 
out side unit radius circle and we will evaluate the residue at Z = Z 2 as 






hence path operator T uu will be 


2 m {Residue at Z = 0 ^Residue at Z =Z 2 } 


Path operator for effective CPA medium will be 


if we write C=Ccp , then 


For CPA medium (+) sign may be (-) sign depending upon the modulus of 
Zi and Z 2 , then 




If we define new variable g 00 such that 


Notice that the path operator T is Hergltoz and g uu is anti Hergltoz. 

The DOS for S-phase shift semicircular model has been calculated by 
Lloyd’s formula 1 as follows: As we are taking binary substitutional 
disordered alloy A X B Y , where X is the concentration of A type atoms in 
alloy and Y is the concentration of B type atoms in alloy. By using equation 
(4.1) the cotangent of phase shift for A and B atoms are 


By using equations (4.5) we can write 


The change in density of states 




where p 0 (E) = free electron density of states per atom, and 

p C p(E) = coherent potential density of states per atom. 
The change in density of states by Lloyd’s formula' is given as 


In above equation, the first term 


where 


In equation (4.29) 


therefore, 


mm 

|||g 



d K 1 1 _ i 


= c B + 


dE 


r« 2.x 2 .x 


U a -U h =x{~~~ ]+r— [ C n-C A ] 

* a 2 .x 


(4.30) 


u = 4-(f l ) = — \XC A + YC B ] 
dE dE L A aJ 


Z 


X Y 

Hh ■ 


r BJ 


2-% 


lxc A + rc B ]+ 


2-X 




U=t 


X Y 
+ ■ 


r, r B 


C i 

■ + • 


2-X 2.x 


(4.31) 


r l -r'=- x [c-c 


(4.32) 


‘a - { H =-x[C A- C H ] 


(4.33) 


By using the equation (4.30), (4.31), (4.32) and (4.33) the explicit form of 


U^i 


IS 


u eS =z 


X Y 

+ ■ 


r A r 8 


+ £ _C) +J£ -C I 


i i 


r r 

1 A 1 B 


(4.34) 


2-X (C A -C B ) 

The change in density of states in equation (4.27) can be written in two parts 


as 


. / 
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tag [First part + Second part] 


where 


First part 


By using equation (4.2 1 ) and (4.34), the first part can be written as 


First part 


Second part 


(Z.C-Lz + B,) dE 


where 


1 


(4.38) 



' ’ X-C-i.x + B. dE K 

With the help of equation (4.2) we can get the energy derivative of B- as 


By using above equation we can write 


Recall the equation (4.37) 


Second part 


Z 2 ,/'(C+B f ) 


2 ./ J '(C + fl,) 


2.X 1 (C + B,,) N 


■b.8(b-B P ).db i 
2 .%\C + b) + N 


2 -x\C + b) 


f°(b)b.db 


■f°(b)db 


( C + b ) 
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Now we can write second part as 
Second part = First integral + Second integral 

where. First integral = — i- f 'L SQML 


We choose the same distribution function f° (b) from equation (4.9) then first 
integral will be 


Putting b = Cos 8 , then 


■ Sin 2 8. Cos 6. d 6 


First integral 


•Sin 2 6. Cos 9. d9 


ince 


First integral 


•Sin 2 9.Cos9.d9 


where 


Using complex analysis it can be written as 


= {' j/(Z)clZ 
where 

(Z 2 -1) 2 (Z 2 +1) 

7 ' Z 3 (Z 2 +2CZ + 1) 

The poles of integrands are Z = 0 of order 3. Other two poles are Z = Zi and 

Z = Z 2 which are roots of Z 2 + 2CZ + 1 

Now, 


(a) Residue at Z = 0 of order 3 

(Z 2 -1) 2 (Z 2 +1) 
(Z 2 +2CZ + 1) 

= 4C 2 -2 


= Limit , 


1 </ 2 


Z-+0 


2! dZ 



(b) Residue at Z = Z 2 


(Z-Z 2 )(Z 2 — l) 2 (Z 2 +1) 
1 Z 3 (Z-Z 2 )(Z -Z,) 



= (Z 2 -Z,)(Zj +Z 2 ) 
= -2(C 2 -l) I/2 (-2C) 


= 4C(C 2 -1) 1/2 

/, = ^-2;n (Residue at Z = 0 + Residue at Z = Z 2 ) 
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Second integral = — — .2[c + (C 2 - 1) 1/2 ] (4.43) 

2-X 


By putting the value of first integral (4.42) and second integral (4.43) in 
equation (4.41), we get the second part as 



55 


Second part 


-1 + 2 C(C ! . 2 [c+(C J - 1 )''=] (4.44) 


The explicit form of ^o(£)can be shown by combining equations (4.36) and 
(4.44) in equation (4.35). 

The change in density of states 


Sp(E) = -lmg 2(C + (C 2 -1) 1/2 )J — + L L (_Ez£LY_L_J_ 

r J 1(^/4 — C B ) r B 


(4.45) 


The change in density of state can be written by using equation (4.23) 


5p{E) = -lmg g 00 — + il| + (i£z£L|f_L_l 

* L ll ^ r J {(c A -c H ) {r A r B 


(4.46) 
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4.3 DOUBLE-PEAK SEMICIRCULAR MODEL 

For this model we will assume (like s-phase-shift semicircular model) 
that the phase shift has resonance behaviour and the cotangent of phase shift 


Ca(bi is as 




where E A( b) is the resonant energy (band center) and Y A{B) is half-width of 

resonance. Another similar approximation is that k-space structure function 
is assumed as 


has no explicit dependence of energy (i.e. independent of energy) where 
X = E U2 and B- is the Fourier transformation of real space structure function 
B nm as given in equation (4.3). 

Now for double peak semicircular model, the explicit form of path 
operator will be similar to s-phase shift semicircular model as 


r 00 = -- f— 8(b-B,)db 
x'N^fC + b K 


(4.47) 


in which the distribution of B f: is 


f°(b) = ±-Y 4 S(b-B R ) 

N ic 


(4.48) 
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and 



Therefore, for double peak semicircular model we choose the distribution 
function f°(b), which must satisfy the above condition. Thus it can be written 


when b > 2 


Then from equation (4.47) , (4.48) and (4.50) 


Let the first term of equation (4.5 1) of T 00 will be as 


First term of T 


Putting (1+b) = X => db = dX 


1 + 1 1-x 2 


fill 
J tt ' r 


x _i ft c+x-i 


Now let C’ = C-1 


■1' , a 

Z*jc C +X 


Now substituting X = Cosd => dX = -SinO.dd 


First term of T Q 


1 y-Co^) L 

X-n J C'+Cos0 


-SinOJO 


1 j Sin 2 6 
X-n ~C+Cos6 


1 2 f Sin 2 6 
2-X-n J C+Cose 


.de { Since f{2n -9) = f{9 ) } 


Now put Z = e‘° => dZ = iZ.dO => d6 


(4.53) 


(Z-J) 


— — f- 

0 V IT J 


2 -^<' c+i z+i 


21 Z 


(Z 2 -1) 

i f L z 2 j 

J (Z 2 + 2CZ'+1) 
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where 


A.x-x " Z 2 (Z 2 + 2CZ'+1) 


The singularities of f(Z) are given by 


Z = 0 (a pole of order 2) 
Z = a (a pole of order 1 ) 
Z = |3 (a pole of order 1) 
we have 


therefore 


Residue of f(Z) at Z = a will be as 


Re 5 f(Z) = Lim(Z-a)f(Z) 
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Lim(Z - a) 


mce a. 


(4.57a) 


Re.v/(Z) 


Similarly, residue of f(Z) at Z=jB will be as 


(4.57b) 


Rev/(Z) 


(ii)Residue at Z = 0 


R tsf{Z) = lim 


4.jjt'Z 2 (Z 2 +2C’Z + 1) 


him — — 
Z-*° A.X-TC 


2(Z - 1).2Z(Z +2CZ + 1)- (Z 2 - 1) .(2Z + 2C 1 ) 
(Z 2 + 2C'Z + 1) 2 ~ ~ 


■i -2 C iC 


A,X-n 1 2xn 


Re.v/(Z) = — 

z.jjt 


Hence by Cauchy Residue theorem 


§ f (Z )dZ = 2n i. [sum of residues inside unit circle] 


2 m 


2-XJt 


-Vc ,2 -i + 


iC 

2-X- k J 


<f/(Z)dZ = -[-C’-Vc : M 


or 


X 

1 


<| f(Z)dz = -[- eWe 2 - i 


if 


if 


(C-1) + V(C-1) 2 -1 


(C-1) + V(C-1) 2 -1 


<1 


>1 


(4.58) 


(4.59) 




Thus, the first term of I 00 will be as 

1 


=> zl °fl.izO±^!L^ = IL cWc , 2 -i] 

^ _ 2 ^r C + 6 


First term of r 


00 


Z 


-(C-1) + V(C-1) 2 -1 


(4.60) 
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where + sign depend on the conditions of equation (4.59) 

Now we take the second term of equation (4.5 1) of 7" 00 will be as 


=> Second term of r°°=-- f- i- - (1 - ^ - db 

Xin C + b 


putting (1-b) =X => db = -dX and let C’ = C+l 


=> 


W c- 


substituting X = Cosd dX - -Sind.dO 


X-7t ; C-CosQ x-x o C'-Cosd 



1 2 f Sin2 ° 
2.x Jc \c-CosB 


.d9 


(4.61) 


putting Z = e' e =>dZ = iZ.dO 


(4.62) 


1 j [(Z-l/Z)/2/f dZ 
2-X-x ' [C'-l / 2.(Z + 1 / Z)] iZ 



Let 


i ( Z 2 -!) 2 

4-X-x Z 2 (Z 2 -2CZ + 1) 


/ (Z 2 -l) 2 

4. Z .*r Z 1 {Z - a){Z - fi) 


(4.64) 


The singularities of f(Z) are given by 
Z 2 |z 2 - 2C" Z + lj = 0 
i.e. Z 2 (Z-a)(Z-/?) = 0 


we have 

cc + /3 = -2C; a.j8 = l & a-/3 = ±2 


a = C'+VC 


/? = C 


* /^« 2 . 


Only one pole among {a,/3) will lie inside the contour of unit circle due to 

conditions a.p- 1. 

Thus the poles lying inside the contour C are 

(i) A simple pole at Z = a or at Z = ft 

(ii) A double pole at Z = 0 
(i) Residue of f(Z) at Z = a 


Re 5 f{Z) = Lim(Z - a) f (Z) 

Z Z 



v i (Z 2 -l) 2 
Lzm(Z a )- 4 ^- Z 2 (Z _ a)(Z _^) 



i (a 2 -!) 2 


4 .x* a*(a-j3) 


{Since a. ft = 1 } 


i ( a i 2 -1/a ) 2 


4 .%Jt {a-fi) 


— —(a - /?) = — - — .2-v/c^-l 
4.J.7T 4./JT 


Re.v f(Z) = — ^ — -V C' 2 — 1 
' 2.J.7T 


or. 


Re.v /(Z) = — — .Vc ,2 -1 

" ' 2.jx 


Z *a 


now. 


ii) Residue at Z = 0 


Res/(Z) = L/m 

z-^-o (2-1)! <3fZ 


1 .-|-[(Z-0) 2 /(Z)] 


i/m— 
? - >0 <afZ 


z 2 . 


2 * 


(Z 2 -!) 2 


4.j./r Z 2 (Z 2 -2C'Z + 1)J 


(4.65) 


l 


Lim 

4.%jc 


2(Z 2 -1).2Z(Z 2 -2C'Z + 1)-(Z 2 -1) 2 .(2Z-2C') 


(Z 2 -2CZ + 1) 2 


i 2C' iC 


4 .xx 1 2.£/r 
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. "• 1 

■ 




M 





hence by Cauchy Residue theorem 

<j> / (Z )dZ =2 jci. [Sum of Residues] 


2m\ 


2-Z-x 


-VcM 


+ • 


iC 


2.x.x 


\f{Z)dZ = --[C'+vCM 


Second term of T 


oo 


\f(Z)dZ = - - |(C + 1) + V(C + 1) 2 - 1 J ; if 


(C+l)+^/(C+l) 2 -1 


<1 


or, 


jf(Z)dZ = - - [(C + 1) - ^(C + 1) 2 - 1 J ; if 


-(C + l) + ^/(C + l) 2 -l 


>1 


(4.67) 


Now using equation (4.60) and (4.67), combining first and Second term of 

rp 00 


,00 


-[(C-l) + iJ(C-l) 2 -1 + (C + 1) + -J(C + 1) 2 - l] 


X 


(4.68) 


If we define new variable g DP such that 


T 00 = 


(4.69) 


X 


Then 


=> g" = |(C - 1) ¥ J(.C - 1) ! - ■ 1 + (C + 1) T -J(C + 1) 2 - 1 1 (4.70) 


Notice that the path operator T 00 is Hergltoz and g DP is anti Hergltoz. 



4.4 DENSITY OF STATES: 

The change in density of states for S-phase shift double-peak semi 
circular model is also calculated by Lloyd’s formula 1 as given in equation 
(4.27), which is recalled: 


Sp(E) = --lmg 

' Tt 


V 00 — 

i ' Jr? ' L 


i ^ 1 


dE / 


;(B«) 


(4.71) 


In above equation let the first term 


T 00 ) = T m .u eff 

\ dE a 


(4.72) 


where 


U* 


X_ Y_ 

L r, + r. 


(i-C) ( C-C ) , 

+ + — — . y 

2-J (Q-Q)^ 1 


J 1_ 

.r, r„ 


(4.73) 


and T 00 for double peak semicircular model is given by equation (4.68) as 


X 


- (C - 1) T J(C - 1) 2 - 1 + (C + 1) + V(C + 1) 2 - 1 J (4.74) 


The change in density of states in equation (4.71) can be written in two parts 
as 


8p{E) = Img [First part + Second part] 

K 


(4.75) 


where. 
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First part 


fc - 1) T -/(C- 1) 2 - 1 + (C + 1) T V(C + 1) 2 - 1 } 



( X y) 

r 

— -f 

l r, rj 


0-c) (C-C) 

+ 1 + — . x\ 

2.x (C A ~C B ) * 


i 

n 

u 


rj 

IJ 


(4.76) 


The second part of equation (4.71) on R.H.S. 


Second part 


l y 1 —(B 1 

N^(t^-B.)dE y q) 


W 


N^i-x-C + Lz-B-) dE q 


_Ly I .—(B.) 

Nfix-C-Lz + Bs) dE K 


_L 

/V 


ZM) 


where, 
/(#*) = 


1 


X.C - i-x + B- dE 


m 


(4.77) 


(4.78) 


With the help of equation (4.2) we can calculate the energy derivative of B - 
as 


dE 


(*,)--* 


5,-s + i 


2-Z 


By using above equation we can write 
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2 -X\C + B R ) 


Now, by using equation (4.77) 


Second part 


2 .% 2 '{C + B r ) 


2 . X 2 \C + B r ) 


2 . Z 2 '(C + B r ) N 


b.5(b- B^-db i 
2 . X \C + b) + N 


2 .z\C + b) 


f\b)db 


(4.80a) 


(C + b) 


We can write second part as 


[ Second integral (I 2 )] (4.80b) 


[First integral (Ii)] + 


Second part 


f{b)b.db 


where, first integral (I j) j 

~oO \ ' / 

We choose the value of f°(b) from equation (4.50) then first integral will be 


b.db + 


69 



Let the first term of above equation (4.81) 


First term of I| 


1 jM +*)*}' 


1/2 


tc_\ C + h 


-.b.db 


putting (l+b) = X=>db = dX 




C + X-l 


'ft- xT JLq-l-f-xSla 

j C + X-l n \ C + X-l 


now let C’ = C -l 


n\ C+X it l C+X 


C+X 


using equation (4.60) 


First term of L 


1 ’fi-X'HXg - [c - 1) + -v/tC-l) 2 -l] (4.82) 

k 4 c+x 1 


Now taking the first part of first term of Ij on R.H.S. of equation (4.82). 


First part of first term of Ii 


L*pt£l 

- i C+X 


K 




Su I i | 
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Now substituting X = CosO => dX = -Sind.de 


Sin9.dd 


C+CosS 


Sin 2 9.Cos0 


C'+Cos0 


Sin 2 O.Cos9 


,dd {Since f(2n -6) = f(6) } 


C'+CosO 


Now putting Z =e‘° =>dZ = iZ.d6 => d6 


i j (Z 2 -l) 2 .(Z 2 +l) 
8/r^Z 3 (Z 2 +2CZ + 1) 


where, 


‘ Z 3 (Z 2 -t-2C’Z + l) 

(Z 2 -1) 2 .(Z 2 +1) 




The poles of integrands are Z = 0 of order 3, other two poles are Z 
the roots of Z 2 + 2CZ + 1. 


= {3, which 
(a) Residue at Z = 0 


on solving we get 


(b) Residue at Z = a 


Re .y f (Z) = Lim(Z - a)f (Z) 

Z»a ‘ 


(g 2 -l) 2 .(a 2 +l) 
a 2 (a 2 -1) 


(a 2 -\).(a 2 +\) 


a 


= (a-l/a) 2 .(a: + l/a) 


(a - fi).{a + p) 


Using equation (4.56) 


(4.86) 



Resf(Z) = (- 2-v/C’ 2 — 1)(— 2C,) = 4 C'VcM 

/f ' 1 

or. 

Re 5' /(Z) = -4 C'VcM 
Using Cauchy Residue theorem 


(4.87) 


First part of first term of Ii= <j> / ( Z ) JZ = 2 % i . [sum of residues] 

= 2^[4C' 2 -2 + 4C’Vc 2 -iJ= 4ot|2C ,2 -1 + 2C , Vc 2 -TJ (4.88) 

From equation (4.82), (4.83) and (4.88) we have 


First part of first term of Ij = — f 

i rr « 


1 *‘ c (l-X 2 ) 1/2 .X 


C'+X 


.dX = — Am 

%7t 


2C ,2 -14 2C , VcM 


or. 


1 A(\-XT-X 
n J C’+JT 


.dX 


l 

2 L 


|2(C - 1) 2 - 1 + 2(C - l)V(C-l) 2 -1 

(4.89) 



. I-.., ... 1 


Therefore, the first part of first integral (fi), of equation (4.81) 


1 ° f (l-(l + 6) 2 ) 1/2 . .. 1 


it 


fiL J±l±LL-.b.db = -42 (C - l) 2 - 1 + 2(C - l)V(C-l) 2 -l]- [(C - 1 ) + y/(C-l?-l 

f 2 C + a l 

(4.90) 


Now, second part of first integral (fi), of equation (4.81). 



Second part of f = - ji 1 (1 h.db 

it i C + b 


putting 1 -b = -X => db = dX 


i T { i -^ 2 }' /2 




n C + l + X 


-.(1 + X).dX 




C + l + X 


,X.dX + - P 

7C ■ jj j C + 1 4 * X 


using equation (4.67) 


Second part of I, = - ft ^— .X.dX + {c + 1) + J(C + 1 ) 2 -l} (4.9 1 ) 

71 4 (C + 1) 4* X 


Now using equation (4.89), second part of first integral (Ij). 
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1 + ril - (} h ) h rih — — — 2 (C + 1) 2 -1 + 2(C + 1)^(C + 1) 2 -1 + j:c + 1) + J(C + 1) 2 -1 

(4.92) 


+ 2 .. 

2T „ C + * 2 


Using equation (4.90) & (4.92), the first integral of equation (4.80) will be as 


I, = *\fM'b.db = -~\l(C + 1) 2 - 1 + 2(C + 1 )V(C + 1) 2 - 1]+ |c + 1) + V(C + 1) 2 - 1 } 
J C + & ? 


-2 


2(C - 1) 2 - 1 + 2(C - l)^(C-l) 2 -! - (C - 1) + V(C-1) 2 -1 


After simplification we have, 


+2 r0 


\i^+b' b ' db = 1 ~ 2C ' 2 ” 4^J- C^C-l) 2 -! 


(4.93) 


The second integral (U) is already solved in equation (4.51) and given in 
equation (4.69) and (4.70) as: 


(c - i)t {J(c-i) 2 -i}+ (c + 1)+ jv/(C+l) 2 -lf 


Therefore, after combining Ii and I 2 in equation (4.80b) we get 
The second part of equation (4.7 1 ) = \^r^- b - db + ^c~+b' db 


+ - 


l —\.c 2 - 1 + C.{](C + \) 2 -l}+ C.|\/(C-1) 2 - l| 

X 

(C-l) + [/(C-l) J -l}+(C + l) + {,/(C + l) J -ll (4.94) 
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Therefore, the change in density of states from equation (4.71) after 
substituting the First part and Second part from equation (4.76) and (4.94) 
respectively. 


8p(E) 


+ 


— Img 
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' c-c 


-((C - 1 ) + ^(C - 1 ) 2 - 1 + (C + 1 ) + V(C + 1) 2 - 1 1 * \x 


( X_ Y_ A 

\T A + r BJ 


-f- 


/ 


\C A -c B/ 


-z- 


r \_ o 

vr«“ t bJ 


2-X‘ 


- { 2 C 2 - 1 + C.-^C + l) 2 -! + C.V(C- l) 2 -l}+ 


l 


2-Z‘ 


|c-l) + V(C-l) ! -l+(C + l) + V(C + l) 2 -l} 


We have defined g'"‘ = ((C-1) + V( C-l) 2 -' + (C + 1) + V(C + 1) J -1 ) in equation 
(4.70), therefore. 


Sp(E) = --lmg\ 

7t 




( X_ 

kTa+Tbj 


-Xg^H^H^Xg") 
£ 2.x 2.x 


1 (gW)f C £ 

* \C,-C 




BJ 


_L_i. 

vr, r, y 
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On solving above equation after neglecting the real part (1/2E), we get the 
change in density of states for double peak semicircular model as 
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From equation (2.23) for S-phase shift double peak semicircular 


where, 


X [X.C A -i.X + Y.C B -i.Y\+ 

[- Z-C + i-X + X-C A ~ i-z}* T 00 * [- X-C + i-X + X-Q 


Using equation (4.69) we have T' 


Then the CPA equation for this model will be as 


=> c = [x.c, + y.c„]+[c-c,].g".[c-c,] (4.98) 

where, 

g n '' = 1 (C - 1) T V(C - 1) 2 - 1 + (C + 1) T ^(C + 1) 2 - 1 


here + sign depends on the conditions. 

Therefore, the change in density of states by Lloyds formula will be as 


Sp(E) = - Im#|c - 1) + #’-l) 2 -i + (C + l) + ^(C + l) 2 -l}* 
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Now for programming we prepare the algorithm first, 


STEP 1 : We declare the cotangent of phase shift for A and B type of atoms 
(i.e. C A & Cb respectively), energy (E), resonance half width for A 
and B type of atom (i.e. T A and r B respectively), band centres of A 
and B (i.e. E A & E B ) and change in density of states (5DOS). 

STEP 2: Input concentration of A (i.e. X), energy (E), resonance half width 
(i.e. r A and f B ) and band centres (i.e. E A & E B ). 


STEP 3: Calculate C A ,C B andC from following formulas 





STEP 5: Calculate g 1 


STEP 6: Calculate [ (C-C)-[C-Cj*g 
derivative (i.e. FDCO). 


STEP 7: Calculate effective (C) by Newton-Raphson’s method 


FDCO 


STEP 8: Repeat the step 3 to 7 till convergence; then calculate the change in 
density of states (i.e. 8DOS) and store energy (E) Vs change in 
DOS (SDOS). 


STEP 9: Increase the energy and repeat step 2 to step 8 for same 
concentration, band centres and resonance half width but for 
different energies. 
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C Programme to calculate the DOS 

COMPLEX CEFF,CA,CB,CBAR,C,FCO,FDCO,GOO,EXPR1 ,EXPR2 
COMPLEX G001,G002,A1,A2,B1,B2 
OPEN (UNIT=3 0,FILE='OUT 1 .CDR') 

CEFF=CMPLX(0.0,0.0) 

X=0.8 

Y=1.0-X 

E=-4.03+CMPLX(0.0,0.00) 

NE=0 
SUM=0.0 
DE=0.03 
50 E=E+DE 
NE=NE+1 
EA=-2.0 
EB=2.0 
GA=1.0 
GB=2.0 

C A=(E A-E)/ GA 
CB=(EB-E)/GB 
CB AR=X* C A+ Y* CB 
C=CBAR+CMPLX(0.0, 0.000000000 1 ) 

IF (DOS .GE. 0.00000001) C=CEFF 

DO 100 1=1,1000 

A1=C-1.0 

B1=CSQRT(C*C-2.0*C) 

A2C+1.0 

B2=CSQRT(C*C+2.0*C) 


IF (CABS(-(C-1 .0)+CSQRT(C*C-2.0*C)) .GT. 1.0) G001=A1+B1 

IF (CABS(-(C-1.0)+CSQRT(C*C-2.0*C)) .LT. 1.0) B1=-B1 

IF (C AB S(-(C+ 1 . 0)+C SQRT (C * C+2 . 0 * C)) .GT. 1.0) G002=A2-B2 

IF (CABS(-(C+1.0)+CSQRT(C*C+2.0*C)) .LT. 1.0) B2=-B2 

G001=A1+B1 

G002=A2+B2 

G00=G001+G002 

FC0=C-CBAR-((C-CA)*(C-CB)*G00) 

FDC0=1 .0+(CA+CB-2.0*C)*G00 
1+(C*(CA+CB))*(2.0+(A1/B1)+(A2/B2)) 
1-(C*C)*(2.0+(A1/B1)+(A2/B2)) 

1-(CA*CB)*(2.0+(A1/B1)+(A2/B2)) 

CEFF=C-(F C0/FDC0) 

IF (CABS(CEFF-C) .LT. 0.0001) GO TO 101 
C=CEFF 

100 CONTINUE 

101 EXPR1=(X/GA)+(Y/GB)+((C-CBAR)/(CA-CB))*((GB-GA)/(GA*GB)) 

EXPR2=EXPR1 *G00 
DOS=AIMAG(EXPR2)/3 .14159 
DOS=ABS(DOS) 

SUM=SUM+DOS * DE 

WRITE (30,*) 'E, DOS, I, SUM', E, DOS, I, SUM 

IF (REAL(E) .LT. 6.0) GO TO 50 

STOP 

END 


. 
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4.6 RESULTS AND DISCUSSION: 

We have calculated the change in density of states ( SDOS ) with 
respect to free electron density of states for A X B Y alloys of resonance half 
widths r., = 1, r# = 1 and 2 , band centers at E A = -2 and E B = 2 for different 
concentrations. 

Figure 4.1 shows the change in DOS of pure B constituent (i.e. X=0.0 
concentration) for unit half width T H = 1 and band centre E B = 2. We see that 
double peaks are forming according to the model at the position of E B ±r H . 
Similarly Figure 4.9 shows the change in DOS of pure A (i.e. X=1 
concentration) for unit half width = 1 and band centre E A = -2. Again we 
see that double peaks are forming according to the model. 


Figure 4.2 shows the change in DOS of alloy for constituent metals 
band centres E A = -2 and E B = 2, resonance half width r„ = r fl = 1 and low 
concentration X=0.1. The double peaks in SDOS are observed around 1 and 
3, which are at the combinations of band centre and half width of B i.e. 
(Eg ± r„). Thus these peaks are arising due to majority of B types of atoms in 
alloy. One single small peak is observed near E A = -2, which is near the band 
centre of A. This peak is arising due to energy levels of low concentrations 
of A types of atoms in alloy. Thus we found that all structures in SDOS are 
arising near band centre and combination of band centres and half widths of 
constituent metals. 
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Figure 4.3 shows change in DOS for same band centres and same 
resonance half width as in figure 4.2 but different concentration X=0.2. We 
see that bigger peaks (double peaks) are near E=1 and 3, which are 
combination of band centres and half width of B. Thus these are arising due 
to majority of B types of atoms in the alloy. One single little flat small peak 
is around E= -2, which arises due to energy levels of minor (20%) 
concentration of atoms in alloy. Again we found that all structures in SDOS 
are arising near band centres and combination of band centres and half width 
of constituent metals. 

Figure 4.4 shows change in DOS for same band centres and same 
resonance half width as in figure 4.2 but different concentration X=0.3. We 
see that major band shows double peaks behavior at E=1 and 3, due to 70% 
majority of B types of atoms in alloy. For this concentration, single peak 
shows more flattering nature due to 30% concentration of A type of atoms in 
alloy and indicate that it will change in double peak after more increase in 
concentration. 


Figure 4.5 shows change in DOS for same band centres and same 
resonance half width as in figure 4.2 but concentration X=0.35. We see that 
majority band show double peak behavior at E=1 and 3, due to 65% majority 
concentration of B types of atoms in the alloy. Minority band shows two 
shoulders near E= -3 and -1, which are combinations of band centres and 
resonance half width of A type of atoms (i.e. E A ± r„ ), thus we found that all 
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structures in SDOS are near combination of band centres and resonance half 
width. 

Figure 4.6 shows change in DOS for same band centres and same 
resonance half width as in figure 4.2 but X=0.4. Now we see that both 
(majority and minority) band shows double peak behavior. The structures in 
minority band are near E = -3 and -1, which are combination of band centres 
and resonance half width of A metal (E A +FJ. Similarly the structures in 
majority band are near combination of band centres and resonance half 
width of B metal. Thus we can say that structures of double peak behavior in 
minority band is clear for more than X=0.35 (or 35% of A type of atoms in 
alloy). 

Figure 4.7 shows change in DOS for same band centres and same 
resonance half width as in figure 4.2 but for high concentrations X=0.5 
(Y=0.5). In first band, we see that the double peak at energies E= -3 and -1, 
which are combination of band centres and resonance half width of A type 
metal (E A ±rj. Similarly in second band we see that the structures at E=1 
and 3 , which are combination {E B ±T B ). Thus, we found that all structures 
in 8DOS are arising near combination of band centres and resonance half 
width of constituent metals. 
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Figure 4.8 shows change in DOS for same band centres and resonance 
half width as in figure 4.2 but X=0.9 (i.e. low concentration 10% of B type 
of atoms in alloy). We see that majority band shows structure near E= -3 and 
-1, which are due to A type of atoms in alloy while minority band shows 
single small peak near E=2, which is due to energy levels of B type of atoms 
in alloy. 


Figures 4.10 to 4.15 show the change in DOS for E A =-2, Eb=2, 
r* = 1 (same as in figure 2) but for double resonance half width r B = 2 and 
different concentrations X=0.0, 0.1, 0.4, 0.5, 0.9 and 1.0. 


Figure 4.10 shows the change in DOS of pure B metals (i.e. X=0.0) 
for double half width r g = 2 and band centres E B =2, we see that double peak 
are at E=0 and 4, which are combination of band centres and resonance half 
width {E b ±r g ). Similarly Figure 4.15 shows the change in DOS of pure A 
metals (i.e. X=l) for unit half width r* and band centres E A = -2. Again we 
see that double peaks are forming at E=-3 and -1, which are combination of 
band centres and resonance half width of A metal {E A ±T A ). 


Figure 4.11 shows the change in DOS for band centres E A =-2 and 
E B =2 resonance half width 1^=1, r B = 2 and low concentration (X=0.10). 
We see that majority band shows the double peaks at E=0 and 4, which are 
combination of band centres and resonance half width of B. Thus these 
peaks are arising due to majority of B type of atoms in alloy. We see one 


single peak between E= -3 and -2, which is arising due to energy levels of 
minority of A type of atoms. 


Figure 4.12 shows the change in DOS of alloy for same parameters as 
in figure 4.11 but for different concentration X=0.4. We see that majority 
band shows structures at combination of band centres and half width of B 
type of atoms in alloy. Minority band shows two structures near E= -3 and 
-1, which are again combination of band centres and half width of A types 
of atoms. 


Figure 4.13 shows the change in DOS of alloy for same parameters as 
in figure 4.1 1 but for high concentration X=Y=0.5. We see that both band 


show clear double peak behavior at the combination of band centres and 
resonance half width of constituent atoms. 

Figure 4.14 shows the change in DOS of alloy for same parameters as 
in figure 4.11 but for different concentration X=0.9. Now we see that first 
band shows peaks at E A = -3 and -1, which are combination of band centres 
and resonance half width of A (E A ±r A ). Another small band shows small 
flat peak between E=1 and 4, which is coming due to energy levels of B type 
of atoms near its band centres Eb=2. 
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Figure 4.9, Change in density of state for rA=1.0,rB=1.0,EA : 

Eb=2.0 for concentration X=1.0 
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Figure 4.11, Change in density of state for rA=1.0,rB=2.0,EA=-2.0 and 

Eb=2.0 for concentration X=0.10 
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Fifure 4.12, Change in density of state for rA=1.0,rB=2.0,EA=-2.0 and 

Eb=2.0 for concentration X=0.40 
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4.7 CONCLUSIONS: 


We have applied KKR-CPA formulation to S-phase shift double-peak 
semicircular model. We have also derived the explicit analytical form of 
change in density of states (5DOS) for double peak semicircular model 
without using step function as in earlier work Yadav 32 , in which it was 
assumed that energy derivatives of step functions was zero. It has been 
found that structures in 8DOS are depending on the combination of band 
centres and bandwidth of constituent atoms. This model can be easily 
applied in the estimation of density of states of alloys in very beginning 
without much computational effort. 




In this thesis, we have discussed and described the applications of 
charge self-consistent KKR-CPA method for calculation of electronic 
structure of disordered alloys using muffin-tin potentials. The charge self- 
consistent KKR-CPA was applied to Cr-Al alloy, which is important bcc 
alloy system. We have calculated self-consistent charge densities, 
component density of states and density of states of Cr x Al Y alloy for 
concentrations X=0.95, 0.85 and 0.75. It is found that charge-self consistent 
KKR-CPA density of states for different concentrations of alloy is in very 
good agreement with reflected X-ray photo spectra (XPS) results. It is found 
that half width of DOS band of these alloys increases approximate 5% at 85 
atomic percentages of Cr with respect to 95 and 75 percent of Cr, which 
explains the change in shape and half width of experimental XPS results. It 
is also found that Cr-based rigid band model fails in the calculation of 
electronic structure of Cr-Al alloys. 

We have applied KKR-CPA formulation to s-phase shift double peak 
semicircular model. We have also derived the explicit analytical form of 
change of DOS (SDOS) for double peak semicircular model. It has been 
found that structures in 5DOS depending on the combination of band centres 
and bandwidth of the constituent atoms. This model can be applied in the 
estimation of density of states of alloy in very beginning without much 
computational effort. 

The KKR-CPA calculations can be improved by going beyond the 
single-site approximation and using the KKR-CPA, which has been 
formulated by Mookergee et al , Gray et al * Razee et al and Rajput et 


al 36 . In this formulation they have used the augmented space formalism, 
which guarantees the herglotz properties of the averaged Green’s function. 

The KKR-CPA calculation can be improved in many other ways. One 
should improve upon the local density approximation, which we had used in 
our calculation. At present many researchers are engaged in this direction 
and have suggested various approximation such as gradient correction, GW 
approximation, weighted density (WD) approximation and self-interaction 

'I'j 

corrected (SIC) approximation etc . We feel that such approximations 
should be first tried and thoroughly tasted on simple systems and pure 

•>o 

metals before these could be applied to disordered alloys. Stauton et al has 
formulated the relativistic version of KKR-CPA theory and can be used to 
improve our KKR-CPA results. However, such corrections in case of Cr-Al 
alloys are supposed to be quite small. The other direction in which the 
theory for disordered alloys can be improved is by going beyond the muffin- 
tin approximation. The KKR-CPA band theory for pure metals has been 
generalized to non-muffin-tin potentials by Brown and Ciftan 39 , Gonis 40 and 
Faulkner 41 . However, their methods have convergence problems as 
discussed by Gonis and Faulkner. To our knowledge the KKR-CPA theory 
has not been formulated for non-muffin-tin potentials so far and is still an 
open problem. 
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